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Quadrupole and octupole softness in the even-even N = Z nucleus 64 Ge is studied on the spherical 
shell model basis. We carry out the shell model calculation using the pairing plus quadrupole ( QQ ) 
plus octupole (OO) interaction with monopole corrections. It is shown that 64 Ge is an unstable 
nucleus with respect to both the quadrupole and octupole deformations, which is consistent with 
the previous discussions predicting the 7 softness and octupole instability. It is demonstrated that 
proton-neutron part Q P Q„ of the QQ interaction is important for the 7 softness or triaxiality. 


PACS: 21.60.Cs, 21.60.Ev, 23.20.-g, 27.50.+e 

Heavy N = Z nuclei with A = 56 ~ 80 show strong 
shape variations such as prolate shape, oblate shape, 
prolate-oblate shape coexistence, and 7 softness, depend¬ 
ing on the mass number. These nuclei lie in transitional 
regions from spherical shape (e.g., 56 Ni Q) to strong pro¬ 
late deformation (e.g., 80 Zr Q). The N = Z = 32 nu¬ 
cleus {gGe 32 is known to be a typical example showing 
7-soft structure in N = Z proton-rich unstable nuclei, 
according to theoretical calculations based on the mean- 
field approximation ||. The calculations predict proba¬ 
ble 7 instability in the ground state, and triaxiality in the 
excited states, i.e., the quadrupole deformation fi 2 ~ 0.22 
and 7 ~ 27°. 

Deformed shell model calculations predict that the nu¬ 
cleon numbers 34, 56, 88, and 134 are strongly octupole 
driving in nuclei where the Fermi surface lies near single¬ 
particle levels with A l = A j = 3 J|]. We can expect that 
nuclei with N = Z near the octupole magic numbers 
exhibit an especially strong octupole effect, because neu¬ 
trons and protons contribute cooperatively. This does 
not necessarily mean a permanent octupole deformation. 
The level pattern of negative-parity states in 64 Ge is not 
a rotational one ||, and the sequence of the 3 _ , 5 _ , and 
7“ levels is irregularly spaced. Thus, 64 Ge is an N = Z 
proton-rich unstable nucleus manifesting a soft structure 
with respect to quadrupole and octupole deformations, 
from experimental and theoretical evidences. In fact, in¬ 
clusion of the 7 deformation improves the E 2 transitions 
of negative-parity states in 68 Ge ||. 

The spherical shell model approach could be more ap¬ 
propriate for describing various aspects of nuclear struc¬ 
ture. It is desirable to add the g 9 / 2 orbital to the full 
Pf shell (/r/2, f>3/2, /5/2, P1/2) for studying both positive 
and negative parity states of 64 Ge, but the shell model 
calculation in this space is impractical at present because 
of the huge dimension. So we restrict the model space 
to the P3/2, f5/2,Pi/2, and g 9 / 2 orbitals, and carry out 
the shell model calculation with the recently developed 
shell model code (b). There are few effective shell-model 
interactions in this model space. 


Recently, an extended P + QQ force |J was applied 
to the /7/2-shell nuclei. This interaction is schematic 
but works remarkably well. The conventional P + QQ 
force was first suggested by Bohr and Mottelson, and 
widely used by Kisslinger and Sorensen, Baranger and 
Kumar, and many authors. Unlike the original applica¬ 
tion to heavy nuclei, the extended P + QQ interaction is 
isospin-invariant In this Rapid Communication, we 
introduce the octupole-octupole (OO) force into the ex¬ 
tended P + QQ force model to describe negative-parity 
states. This interaction is quite useful for studying not 
only the 7 softness but also the octupole instability men¬ 
tioned above. The QQ and OO forces are the long-range 
and the deformation-driving part of the effective interac¬ 
tion. Contrary to this, the monopole pairing force can 
be associated with short-range force, and restores the 
spherical shape. Thus, the competitions among the QQ, 
OO, and monopole pairing forces are expected to be im¬ 
portant for shape transitions of quadrupole and octupole 
deformations in 64 Ge. The P + QQ force with the OO in¬ 
teraction will be suitable for studying the monopole pair¬ 
ing, quadrupole, and octupole correlations. Recently, the 
fp shell model calculation 1101 with the FPD6 interaction 
jlTf has been performed in "^Ge as a test case for quan¬ 
tum Monte Carlo diagonalization (QMCD) method. The 
projected shell model calculation |lj] has been performed 
in 64 Ge modifying the standard Nilsson parameters. 

Since protons and neutrons in the N = Z nuclei oc¬ 
cupy the same levels, one would expect strong proton- 
neutron ( p-n ) interactions [ p~5| . In particular, one of the 
most interesting questions in the study of nuclear struc¬ 
ture is what roles the p-n interaction play in the nuclear 
deformation. The long-range p-n isoscalar (T = 0) in¬ 
teraction between valence nucleons has been suggested 
to be a source of the nuclear deformation jl(|. On the 
other hand, the isoscalar QQ interaction used in the 
P + QQ force model has very strong p-n component 
QpQn , which gives rise to nuclear quadrupole deforma¬ 
tion [|l7|]. The QpQn interaction is expected to be impor¬ 
tant for quadrupole collectivity in 64 Ge. The rotational 
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behavior of T = 0 and T = 1 bands in the odd-odd 
N = Z nucleus 62 Ga is recently studied using the spheri¬ 
cal shell model and the cranked Nilsson-Strutinsky model 



In order to study the octupole correlation, let us in¬ 
troduce an isoscalar octupole interaction Hoo with the 
force strength y 3 to the extended P + QQ model || with 
monopole corrections f/ 1 ™ rr : 

H = H 0 + H Po + Hp 2 + H qq + H 00 + H c ™ 

= ^ gj PjmIk,P jm1k 

a J—0,2 Mk, 

E « : QLQ2M : £ X 3 : 0\ m Oom : +H%“, 

M M 

(i) 

where e Q is a single-particle energy, Pjmtk. is the pair 
operator with angular momentum J and isospin T, and 
<? 2 M (0 3 m) is the isoscalar quadrupole (octupole) op¬ 
erator. Due to the isospin-invariance, each term of the 
above Hamiltonian includes p-n components, which play 
important roles in N — Z nuclei. 

We carried out shell model calculations in a model 
space restricted to the 2 p 3 / 2 , 1/5/2, 2pi/2, and lgg / 2 or¬ 
bitals (called pfg- shell henceforth). The model assumes 
a closed IgN^s core and does not allow for core breaking. 
The neutron single-particle energies of 2 p 3 / 2 , I/5/2, ^Pi/2, 
and l^g /2 in this pfg- shell region can be read from the 
low-lying states of 57 Ni, because the low-lying states of 
57 Ni are well characterized as pure single-particle levels 
when 56 Ni is a closed shell core. The adopted single¬ 
particle energies relative to the 2 p 3 / 2 are e p3 / 2 = 0 . 0 , 
££ 5/2 = 0.77, £pi / 2 = 1.11, and e g g /2 = 3.70 in MeV 
[|19| . Since the above Hamiltonian is assumed to be an 
isospin-invariant, the proton single-particle energies are 
taken as the same values as the neutron single-particle 
energies. The force strengths of the extended P + QQ 
interaction are taken so as to reproduce the energy levels 
of low-lying states in 64 Ge as follows: 

go = 0.426(42/H), g 2 = 0.274(42 /A) 5 / 3 , 

X 2 = X 2 (42 /A) 5/3 /b 2 = 0.567(42 /A) 5 / 3 /b 2 , ( 2 ) 

X3 = x°(42 /A) 2 /b 6 = 0.275(42 /A) 2 /b e , 

where go, g 2 , X2> and X 3 are the monopole pairing, 
quadrupole-pairing, QQ, and OO force strengths, respec¬ 
tively. We adopt the harmonic-oscillator range parameter 
b ~ H -1 / 3 , the effective charge e p = 1.50e for proton and 
e n = 0.50e for neutron. We adjust phenomenologically 
force strengths of several monopole corrections so as to 
approximately reproduce the low-lying energy levels of 
64 Ge. These force strengths can also reproduce the low- 
lying energy levels of 58 “ 66 Ni, 60 " 64 Zn, 66 Ge, and 68 Se. 

In Fig. 1, calculated energy spectra are compared with 
experimental data for 64 Ge. Two side bands are shown 


in addition to the ground-state band, i.e., positive-parity 
band on the band head 2 + and negative-parity band on 
the band head 3 _ . The calculations reproduce the ob¬ 
served three bands at good energies. The agreement be¬ 
tween theory and experiment for the ground-state band 
up to spin I = 8 is good. The calculated B(E 2) value be¬ 
tween the ground state and the first excited I = 2 + state 
is B{E2\2 J [ —> 0j l ")=245.3 (e 2 fm 4 ) corresponding to the 
quadrupole deformation (3 ~0.2. This value is consistent 
with the predictions of f3 ~0.22 by Moller and Nix |2C| ] 
and of (3 ^0.22 by Ennis et al. ||, and is comparable to 
the experimental data 12 W.u. of 66 Ge nucleus. 
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FIG. 1. Comparison of experimental and calculated energy 
levels of 64 Ge. The arrows designate E2 transitions with the 
calculated B(E2) values indicated by their widths. 

The calculated occupation numbers of the p 3 / 2 , f 3 /2, 
Px/ 2 , and gg / 2 orbitals in the ground-state band are 3.6, 
3.2, 0.6, and 0.6 on the average, respectively. More than 
four nucleons are excited from the unperturbed configu¬ 
ration (p 3 / 2 ) 8 - The full fp shell model calculation 0 
with the FPD 6 interaction |ll[ using the QMCD method 
has recently been performed for low-lying I = 0 ^, 2f, 2%, 
and 4+ states of positive parity in 64 Ge. The FPD 6 cal¬ 
culation predicts the deformation /?2 ~ 0.28 which is 
somewhat larger than the other predictions (3 2 ~ 0.22, 
and gives the triaxiality 7 ~ 27° which is consistent with 
the others predictions. The FPD 6 interaction seems too 
strong to yield appropriate collectivity [^ 2 |, due to its 
drawback [|l3|. This can be seen from the occupation 
numbers of f 7 / 2 ,P 3 / 2 , /s/ 2 , and pi /2 which are 15.1, 2.6, 
5.5, and 0.8, respectively. Two more nucleons are jump¬ 
ing to the orbitals above p 3 /2 as compared with our re¬ 
sult. The stronger collectivity caused by the FPD 6 inter¬ 
action is probably attributed to the mixture of the three 
orbitals (/V/ 2 , p 3 /2, and f 3 / 2 ) due to the large matrix el¬ 
ements between {/t/ 2 ,P 3 / 2 ) and f 3 / 2 - This is the reason 


2 





why B(E2-,2+ —► 0+) = 5 x 10 2 (e 2 fm 4 ) obtained in Ref. 
Jl0[ is almost twice that of ours 245.3 (e 2 fm 4 ). 
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FIG. 2. The excitation energies of the first and second ex¬ 
cited 2 + states in the upper figure and the 0 + ground-state 
energy in the lower figure, as a function of the quadrupole 
force strength. The arrows indicated by their widths desig¬ 
nate E 2 transitions with the calculated B(E 2) values. 

The ratio of excitation energies for the 4f state and the 
2 ^ state gives additional information with respect to the 
shape. This ratio is 3.33 for a rigid rotor, 2.0 for a pure vi¬ 
brator, and ~2.4 for a 7 -soft nucleus. Therefore, the ratio 
2.65 in the calculation indicates 7 softness. We can also 
see the 7 -soft nature in E2 transitions. The B(E 2) value 
of 2 ^ —► 2f transition is larger than i?(i? 2 ; 2 j —> 0 ^), 
and the ratio B{E2\2\ —> 2\)/B(E2\2\ —> oj) is ~ 27, 
corresponding to 7 ~ 26° in the Davydov model jnj. 
The present result is in agreement with the triaxiality 
estimated from the experimental data and the other the¬ 
oretical models, though it is difficult for the shell model 
calculation to definitely discuss 7 softness and triaxiality. 

On the other hand, the low-lying negative-parity states 
with I = 3 _ ,5 _ ,7~, and 9“ are nicely reproduced. The 
octupole force newly introduced to our model plays an es¬ 
sential role in describing the negative-parity states. The 
value of B(E2; 5“ —> 3 _ ) is small, while the B(E2) values 
between higher states above 3~ are large. This indicates 
a different structure of 3“ from the higher negative-parity 
states. The same is known from the occupation numbers 
°f P3/2> /5/2J Pi/2) and 99/2 orbitals. Their values of the 


I = 3“ state are, respectively, 4.0, 1.2, 0.6, and 2.2, while 
the average values of the negative-parity states above 3“ 
are 3.4, 2.7, 0.5, and 1.4. More nucleons are jumping 
from the unperturbed configuration (7*3/2) 7 (59/2 ^ in the 
3“ state as compared with the states above 3 _ . The 
I = 3“ state is mainly a coherent mixture of neutron 
and proton (7*3/2/5/259/2)3- configurations. The present 
result is consistent with the calculations by Petrovici and 
Faessler j|). In Fig. 1, the calculated 11“ and 13~ lev¬ 
els are considerably higher than the observed excitation 
energies. This may suggest that excitation from the 
orbital in 56 Ni core should be taken into account for these 
states. 

Figure 2 shows the excitation energies of the first and 
second 2 + states and the 0 + ground-state energy as a 
function of the quadrupole force strength \ 2 - The other 
force strengths are fixed to the values of Eq. (2). We 
can see a level crossing of the first excited 2j*~ state and 
the second excited 2 \ state near the crossing strength 
X 2 cr ~0.45 MeV. When the quadrupole force strength is 
apart from X 2 cr; the second excited 1 = 2^ state rises 
abruptly in energy, while the first excited 7 = 2 / state 
does not change largely. This means that 7 softness or 
triaxiality is realized beyond the crossing point. The 
ground-state energy rapidly decreases beyond this cross¬ 
ing point, while it is almost constant for 0 < X 2 < 0-45 
MeV. 
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FIG. 3. Excitation energies of the first and second excited 
2 + states as a function of the p-n quadrupole force strength 
X 2 P n ■ The arrows designate E 2 transitions with the calculated 
B(E2) values indicated by their widths. 

There seems to be a critical point of phase transition 
with respect to the quadrupole correlation. Indeed the 
B(E2;2f — * Oj 1 ") values dramatically increase when the 
quadrupole force strength goes beyond the critical point. 
The force strength adopted in Fig. 1 is X2=0-567 MeV 
denoted by a solid circle in the lower part of Fig. 2, 
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and its position lies in the deformed region but near the 
critical point. We can also see that the B(E2 ; 2% —*> 2+) 
value is large for x§ > 0-45 MeV. This is consistent with 
the 7 softness or the triaxiality discussed above. 
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FIG. 4. Excitation energies of the negative-parity states 
(in the upper figure) and the 0 + ground-state energy (in the 
lower figure) as a function of the octupole force strength. The 
solid circles show the energies for the adopted force strength. 

As mentioned earlier, it is very interesting to study 
roles of the proton-neutron part of the QQ force ( Q P Q n )• 
We make the study by varying the strength of Q p Q n 
(X 2 P n) an< 4 keeping the other force strengths of Eq. (2). 
This is effective in seeing the dependence of quadrupole 
deformation on Q p Q n , though the Hamiltonian without 
the isovector type of QQ force stops beeing isospin invari¬ 
ant. Figure 3 shows the excitation energies of the first 
and second 2 + states and the B(E 2) values as a function 
of X 2 P n- The first I = 2+ state is flat in energy with 
respect to the force strength X% n ■ However, the second 
excited I = 2^ state strongly depends on x% P m and be¬ 
comes lowest near X 2 P n ~ 0-35 MeV. In the isoscalar QQ 
force, P n ^ equal to the magnitude of proton-proton 
(pp) and neutron-neutron (nn) force strengths due to the 
isospin-invariance, i.e., x% n = X% P = xLn- The force 
strength X 2 P n = 0.567 MeV which is adopted in Fig. 1 
leads to the large B(E2) value B(E2\ 2+ —> 0^)=245.3 
(e 1 2 3 fm 4 5 6 7 8 ) corresponding to the quadrupole deformation 
/3 ~0.2, as mentioned above. We can see, however, that 
B(E2\ 2+ —> 0+) is very small for 0 < x% n < 0-35 MeV. 


This suggests that a phase transition occurs near the crit¬ 
ical point X 2 P ncr ~ 0.35 MeV. The B(E2) value of the 
2j —> 2+ transition also becomes large for X% n > 0.35 
MeV, which results in the expected triaxiality 7 ~ 26° 
when X2 pn = 0.567 MeV. Thus, the Q p Q n interaction 
plays an important role for the 7 softness or the triaxial¬ 
ity in 64 Ge. 

Let us lastly study the negative-parity states as a func¬ 
tion of the octupole force strength in Fig. 4. The other 
force strengths are fixed to those of Eq. (2). As men¬ 
tioned above, the I = 3J~ state is very collective, and 
the excitation energy of the I = 3)" state decreases as 
the octupole force strength increases until x° ~0.3 MeV, 
and increases as it goes beyond this point. The other 
negative-parity states have an insignificant dependence 
with respect to the force strength x '3 until the critical 
point, and increase for X 3 > 0.3. The ground-state en¬ 
ergy is almost constant for 0 < x' 3 cr < 0.3 MeV and 
decreases quickly when going beyond this critical point. 
It seems that a phase transition occurs near the criti¬ 
cal force strength x' 3 Cr ~0.3 MeV. The octupole force 
strength Xg=0.275 MeV adopted in Fig. 1 is very close 
to the critical point x' 3 cr- Thus, 64 Ge seems to be near 
an octupole instability. 

In summary, we have studied quadrupole and octupole 
correlations in the even-even N = Z nucleus 64 Ge by 
means of spherical shell model calculations. The P + QQ 
force model including octupole interaction and monopole 
corrections, which is schematic but realistic, is adopted 
for describing the quadrupole and octupole correlations. 
It is shown that 64 Ge is an unstable nucleus with re¬ 
spect to both the quadrupole and octupole deformations. 
The present results reveal that the p-n QQ interaction 
(QpQn) induces an onset of quadrupole deformation and 
7 softness. It can be expected to play an important role 
in the prolate-oblate shape coexistence of the neighboring 
even-even N = Z nucleus 68 Se, which has been recently 
observed p2| , p3[ . 
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